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ABSTRACT 

We present a new procedure to measure the large-scale carbon monoxide (CO) emissions across cosmic 
history. As a tracer of the large-scale structure (LSS) itself, the CO gas content as a function of redshift 
can be quantified by its three-dimensional fluctuation power spectra. Furthermore, cross-correlating 
CO emission with other tracers of LSS offers a way to measure the emission as a function of scale and 
redshift. Here we introduce the formalism and model relevant for such a cross-correlation measurement 
between CO and other LSS tracers, and between different CO rotational lines. As an illustration, we 
propose a novel use of Cosmic Microwave Background (CMB) data and attempt to extract redshiftcd 
CO emissions embedded in the Wilkinson Microwave Anisotropy Probe (WMAP) dataset. We cross- 
correlate the all-sky WMAP7 data with LSS data sets, namely, the photometric quasar sample and 
the luminous red galaxy (LRG) sample from the Sloan Digital Sky Survey (SDSS) Data Release 6 and 
7 respectively. We are not able to detect a cross-correlation signal with either CO(l-O) nor CO(2-l) 
lines at different redshifts, mainly due to the instrumental noise in the WMAP data. However, we 
are able to rule out models that are more than three times greater than our more optimistic model. 
We discuss the cross-correlation signal from the thermal Sunyaev-Zeldovich (tSZ) effect and dust as 
potential contaminants, and quantify their impact for our CO measurements. We discuss forecasts for 
current CMB experiments and a hypothetical future CO focused experiment, for which we propose to 
cross-correlate the CO temperature data with the Hobby-Eberly Telescope Dark Energy Experiment 
(HETDEX) Lya-emittcr sample, for which a signal-to-noise ratio of 58 is possible. 
Keywords: intergalactic medium; cosmology: observations; diffuse radiation; large-scale structure of 
the Universe 



1. INTRODUCTION 

The correlation of redshifted emission line spectra has 
been proposed as a potentially powerf ul probe of large- 
scale structure (LSS) at high redshifts (ISuginohara et al.l 
fl99llRighi etalll200aiVisbal fc Loebll2010D . While tra- 
ditional LSS probes such as galaxies and quasars for the 
foreseeable future will be limited to z < 2 and z < 4, re- 
spectively, mapping emission lines can probe out to z > 
6, potentially mapping the Epoch of Rcionization (EoR). 
In these very high-redshift regions, emission lines can also 
be cross-correlated with the 21-cm line from the inter- 
galactic medium (IGM), constituting a very complemen- 
tary probe to 21-cm alone or extragalactic emission lines 
alone, whose autq - spectra are plagu ed with foregrounds 
(|Gong et all 1201 H iLidz et all 1201 11 ). Three extragalac- 
tic emission lines of interest i n these types of studies 
inclu de CO (IRighi et al.l l2008t ICarillil [20111 IGong et all 
[20111: ILidz et all l201lT l Clf lGong et al.N2012D . and HI 
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In this paper, we will focus on the 



(IGong et al.ll20Tl . 
CO lines. 

CO is readily produced from carbon and oxygen in 
star-forming regions. The CO molecule exhibits several 
rotational transitional states J — >• J — 1 with line fre- 
quencies vj = Juqo and uco = 115 GHz. CO lines have 
been studied both as foreground contaminan ts to Cosmic 



Micro wave Background (CMB) observations ( Righi et all 
2008 ) as well as LSS tomography probes (jGong et all 



2011 



|Lidz_el_alJ|2Qll|) . In this pa per w e will follow the 
formalism presented in lLidz et al.l ()2011l ) (hereafter Lll) 
and extend it to lower redshifts, producing two mod- 
els that make different assumptions. The two CO rota- 
tional lines we seek to constrain are CO(1-0) and CO(2-l) 
because they are typically the brightest and both their 
emitting and redshifted frequencies conveniently fall in 
the frequency range of CMB observations, where vast 
amounts of data and expertise have been cultivated over 
the years. The multiple emission lines of CO due to mul- 
tiple rotational states allow interlopers to be identified 
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more easily, giving CO an advantage over other emis- 
sion lines as cosmological probes. These observational 
prospects have created considerable interest in building 
future dedicated experiments, even if the exact strength 
of the signal at high redshift and on large scales is still 
rather uncertain. As such, we investigate what can be 
probed with current experiments. For the LSS study of 
interest, we work in the regime of "intensity mapping" , 
measuring aggregated CO emissions associated with the 
large-scale structure rather than with individual galax- 
ies, as an effic i ent w ay of extracting the faint signals (e.g., 
IChang et all This approach is complementary 

to high-resolution CO observations enabled by the At- 
acama Large Millimeter/submillimetcr Array (ALMAEI) 
and the Jansky Very Large Array (JVLAEIh for exam- 
ple, where individual galaxies have been mapped to high 
rcdshifts; their small (sub-arcminutc) field-of-view is not 
well suited for large scale surveys. 

One probe that could possibly constrain CO line emis- 
sion is in fa ct the Wilkinson Mi crowave Anisotropy probe 
(WMAP) (jJarosik et all 120111 ). With band frequencies 
ranging from 23-94 GHz, WMAP should contain emis- 
sion up to z ~ 4 from CO(l-O) and up to z ~ 7 
from CO(2-l). Constructing an auto-correlation spec- 
trum in each band is very difficult due to contamina- 
tion from other astrophysical processes that motivated 
WMAP in the first place such as the Cosmic Microwave 
Background (CMB) or galactic emission. However, a 
cross-correlation with a different LSS tracer would not 
have the same foregrounds and thus be immune to these 
contaminants. A suitable choice would be to cross- 
correlate the WMAP bands with photometric quasars 
in its corresponding CO redshift range. Specifically, we 
choose the photometric quasar sample constructed from 
the Sloan Digital Sky S urvey (SPSS') (lYork et al J 12000ft 



D ata Release 6 (DR6) ( Adelman-McCarthv et al 



2008) 



in lRichards et al.1 ()2009l ). This sample contains quasars 
with redshifts as high as z = 6.1, which covers most of 
WMAP's CO redshift range, allowing us to perform this 
cross-correlation for both CO(1-0) and C0(2-l) lines. 
Note that this redshift range could allow us to search 
for higher CO lines. CO(3-2) would be seen in V and W 
bands, while CO(4-3) and CO(5-4) would be seen only in 
the W band. Wc leave searches for these lines to future 
work as they are expected to be dimmer. At redshifts 
z < 1, we can cross-correlate the WMAP W band with 
the SDSS Data Release 7 (DR7) (|Abazaiian et all 12009ft 
spectroscopi c Lum inous Red Galaxies (LRGs) sample in 
iKazin et all (|2010ft . 

In this paper, we propose cross-correlating brightness 
temperature maps with quasars (QSOs) and Luminous 
Red Galaxies (LRGs) to detect fluctuations in CO at 
high redshifts. We first see what can be done with cur- 
rent data by cross-correlating WMAP temperature maps 
in each band with photo-quasars from SDSS DR6 and 
LRGs from SDSS DR7 in the band's CO redshift range. 
We attempt to remove CMB fluctuations from the maps 
to isolate the CO signal (see Sec. 0]for details). We re- 
port no detection of a cross-correlation between quasars 
or LRGs and CO temperature due to extragalactic CO 
line emission, although we place an upper limit that rules 

1 https://almascience.nrao.edu/ 

2 https : / /saf e .nrao . edu/evla/nova/index . shtml 



out any brightness level much greater than our model. It 
appears that the weakness of these constraints is mainly 
due to instrumental noise in the WMAP temperature 
maps, although the density of quasars plays a lesser role. 
We do not suspect quasar photo-z errors to play a role 
since the errors are Az ~ O(0. 1) while the redshift bins 
we use are of order unity. To detect these lines in the 
future, we look towards a hypothetical future CO exper- 
iment to measure CO fluctuations for the line CO(l-O) at 
z = 3 over the range 2.9 < z < 3.1 with 20 redshift bins. 
We forecast that this experiment cross-correlated with 
the full spectroscopic quasar survey fro m the Baryon Os- 
cillation Spectroscopic Survey (BOSS) ([Eisenstein et al.l 
[20rllRoss et~aTJ[2012l) will measure the cross-correlation 
amplitudes with signal-to-noise ratios (SNR) as high as 4 
for each redshift bin and 20 for the entire redshift range, 
depending on the model. The CO experiment cross- 
correlated with the Hobb y-Eberly Telesco pe Dark Energy 
Experiment (HETDEX) ((Hill et al.ll2008h could do signif- 
icantly better with an SNR of 13 for each redshift bin and 
58 for the entire redshift range. We also set constraints 
for current and future ground-based CMB polarimetcrs 
like SPTPol or ACTPol. 

The plan of our paper is as follows: in Sec. [21 we de- 
scribe the data products we us e, including the W MAP 
7-year observations (WMAP7) (jJarosik et al.ll201lD tem- 
perature maps, the photo-quasar maps from SDSS DR6, 
and the LRG maps from SDSS DR7. In Sec. El wc derive 
the form of the CO-LSS cross-power spectra along with 
statistical errors for both models. In Sec. 21 we describe 
the estimator for the cross-correlation, presenting mea- 
surements in Sec. [5] Finally, we present forecasts for a 
possible future CO experiment as well as ground-based 
experiments in Sec. [6] and conclude in Sec. [7] Wherever 
not explicitly mentioned, we assume a flat ACDM cos- 
mology with parameters compatible with WMAP 7. 

2. DATA 

2.1. WMAP Temperature 

We attempt to extract redshifted CO e mission associ- 
ated with LSS embedded in the WMAP7 (jJarosik et al.l 
1201 lf PI temperature maps by means of cross-correlation. 
Primordial CMB and Galactic foregrounds which con- 
tribute to the WMAP temperature maps will not cor- 
relate with LSS; however, CMB secondary anisotropy 
signals, in particular the thermal Suny aev-Zeldovich ef- 
fect (tSZ) IjSunvaev fc Zeldovichl 119801 ) . are expected to 
correlate with LSS. TSZ has be e n det ected using simi- 
lar data sets in IChatterjee e t al.1 (j2010j ) and we will dis- 
cuss it more in Sec. 15.31 Radio and dust emission also 
should correlate with LSS, and we discuss its contribu- 
tion in Sec. l5.3l as well. WMAP has 5 temperature bands: 
K (23 GHz), Ka (33 GHz), Q (41 GHz), V (61 GHz), 
and W (94 GHz) with bandwidths of 5.5, 7.0, 8.3, 14.0, 
and 20.5 GHz, respectively. We do not use the K map 
since it is dominated by Galactic emission; thus we are 
left w ith 4 bands. We use the HEALPix (jGorski et al.l 
120051 ) N res — 9 maps with the KQ85 mask to remove the 
brightest resolved point sources and the bright Galactic 
emission regions near the galactic plane. These cuts leave 
each map containing 2,462,208 N res = 9 pixels, covering 
32,289 deg 2 (78.3% of the sky). 

3 http:/ /lambda. gsfc.nasa.gov/product/map/currcnt 
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One subtlety that must be taken into account is that 
the temperature values in the WMAP data products 
are perturbations of physical temperature, assuming a 
mean temperature equal to the CMB temperature. The 
CO temperature, however, is a brightness temperature. 
These two temperature measurements diverge at the 
higher frequency bands, mainly the V and W bands. 
Therefore, before performing our analysis we must con- 
vert the CMB physical temperature perturbations to 
brightness temperature perturbatio ns according to the 
formula in iGawiser fc Smootl (| 19971 ) 



FWHM Gaussian) n sta rs > 562 stars /deg 2 (twi ce the 

The 



6T b 



x 2 e x 



:5T 



(1) 



where 8Tb is the brightness temperature, ST is the tem- 
perature map from the WMAP data product, and x = 
hv/kTcMB with Tqmb = 2.725A' and v being the band 
frequency. The prcfactor for ST has a different value for 
each WMAP band, ranging from 0.97 for the Ka band to 
0.80 for the W band. Note that this expression neglects 
the band width. 

2.2. SDSS Data 

We use photometric quasars from the SDSS DR6 
(jAdelman-McCarthv et~atl 12008ft a nd LRGs from the 
SPSS DR7 (jAbazajian et al.1 120091 ) to trace the mat- 
ter density and construct its angular power spec- 
trum. The SD SS consists of a 2.5-m telescope 
(IGunn et al.l 120061) with a 5-filter (ugriz) imaging cam- 
era ([Gunn et al.l |1998[ ) and a spectrograph. Auto- 
mated pipelines are resp onsible for the astrometric so- 
lution (IPier et all 120031) and pho t omet ri c calibration 



(iFukugita et al J 119961: iHogg et all I2MH iTucker et all 
120061: iPadmanabhan et all I2008D . Bright galaxies, 
LRGs, an d quasars ar e selec t ed for follow-up spec- 
troscopy (IStrauss et all I2002t lEisenstein et"all 120011: 
IRichards et alll2002l:IBlanton et al.ll2003l) . Thedata used 
here from DR6 and DR7 were acquired between August 
1998 and July 2008. 

2.2.1. SDSS quasars 

We use the photom etric quasar catalog composed by 
IRichards et all (|2009D (hereafter RQCat). The entire 
catalog consists of 1,172,157 objects from 8417 deg 2 
on the sky selected as quasars from the SDSS DR6 
photometric imaging data. Wc limit our dataset in this 
analysis to UV-excess quasars (catalog column uvxts 
= 1) because they have a higher selection efficiency. 
We also require the catalog column good > to reject 
objects that are likely stars. For the survey geometry 
we use the union of the survey runs retrieved from the 
SDSS DR6 CAS server. We omitted runs 2189 and 
2190 because many objects in these runs were cut from 
RQCat. This mask was pixelized usin g the MANGLE 
softw are (Hamilt on k, Tegmarkl 120041 : iSwanson et all 
I2008D . We pixelize the quasars as a number overdensity, 
5 q = {n q — n) /n, onto a HEALPix pixelization of the 
sphere with N res = 9, where n is the pixel's number 
of quasars divided by the pixel's survey coverage w 
and n = w i)- We then reject pixels 

with extinction E{B — V) > 0.05, full widths at half- 
maximum of its point-spread function (PSF) FWHM 
> 2 arcsec, and stellar densities (smoothed with a 2° 



average stellar density) , similar to IHo et all 
extinction cut is very important because a high extinc- 
tion affects the u band, which is crucial to identifying 
quasars. Also, since stars tend to be misidentified as 
quasars, it seems prudent to cut regions with high stellar 
densi ty. We implement these cuts using dust maps 
from iSchlegel et all (fl99l and stars (18.0 < r < 18.5 ) 
from the SDSS DR6 (jAdelman-McCarthv et al.l 120081) . 
We also reject pixels for which the survey region 
covers less than 80% of the pixel area. In addition, 
RQCat contains regions that seem undersampled. 
We excise angular rectangles around these regions 
to remove them from the data. The angular rect- 
angles in celestial coordinates that were removed are 
(a, S) = (122 o -139°,-1.5 o -(-0.5)°), (121°-126°, 0°-4°), 
(119°-128°, 4°-6°), (111°-119 , 6°-25°), 

(111.5°-117.5 , 25°-30°), (110°-116°, 32°-35°), 

(246°-251°, 8.5°-13.5°), (255°-270°, 20°-40°), 
(268°-271 ,46°-49 ), and (232°-240°, 26°-30°). Finally, 
we cut pixels that appeared to have severe photometric 
calibration errors (> 5 mag). After these cuts, the sur- 
vey region comprises 534,564 N res = 9 pixels covering 
7010 deg 2 . 

2.2.2. SDSS LRGs 



We use the LRG catalog composed by iKazin et al.1 
(pOlOl ). LRGs are the most luminous galaxies in the uni- 
verse, making them important for probing large volumes. 
They are also old stellar systems with uniform spectral 
energy distributions and a strong discontinuity at 4000 
A, which enable precise photometric redshifts. The LRG 
catalog consists of 105,623 spectroscopic LRGs from red- 
shifts 0.16 < z < 0.47. We do not make any alterations 
to this catalog. Our survey region comprises approxi- 
mately 638,583 N res = 9 pixels covering 8374 deg 2 . 

2.3. WMAP-QSO Cross-Data Set 

We intersect the pixel sets from the temperature and 
quasar maps to produce two sets of maps that can be 
cross-correlated. This operation leaves each map with 
441,228 N res = 9 pixels covering 5786 deg 2 with ~ 7 
arcmin. pixel resolution. Since each WMAP band probes 
a separate redshift range for each of the CO emission 
lines, we also divide the quasar map into 8 maps, two for 
each of the WMAP bands. Note that some of the bands 
for CO(1-0) will intersect in redshift with other bands for 
CO(2-l). We list the properties of the 4 WMAP maps 
and 8 quasar maps in Table[T]as well as the probed spatial 
scale determined by our pixel size. Note that the Ka(2- 
1) redshift range exceeds the redshifts of the quasars; 
therefore, wc will not determine any limits on the CO (2- 
1) line with the Ka band. 

2.4. WMAP-LRG Cross-Data Set 

We also intersect the areas of the temperature and 
LRG maps for cross-correlation, with 619,708 N res = 9 
pixels covering 8126 deg 2 . Note that the WMAP W band 
is the only band that overlaps with the LRG sample, and 
this is true only for the CO(1-0) line. Thus, we will only 
get one constraint from this cross-correlation. However, 
the number of LRGs in this redshift range is much more 
than the number of quasars, so we expect a more signif- 
icant constraint than that from the quasars. 
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Table 1 

WMAP rcdshift bins for CO emission lines J = 1 — > and J = 2 — > 1 with the quasar (QSO) counts from the QSO map (T x QSO 
intersected map). For reference, we write in parentheses the transverse scale corresponding to the pixel scale in Mpc/h for each z band 

probed. As stated above, the total number of QSOs in DR6 is 1,172,157. 



Band 


z(1^0) QMpc/JiJ) 


AqsoU -> 0) 


z(2 -► 1) ([Mpc/fcJ) 


JV QSO (2 -> 1) 


Ka 


2.151-2.898 (9) 


90,780 (74,395) 


5.302-6.796 (13) 


80 (63) 


Q 


1.547-2.121 (7) 


146,111 (121,297) 


4.094-5.242 (11) 


573 (466) 


V 


0.691-1.130 (5) 


62,172 (51,818) 


2.382-3.260 (9) 


27,032 (22,180) 


w 


0.103-0.373 (1) 


42,184 (34,400) 


1.206-1.746 (6) 


140,819 (116,791) 



3. CROSS-CORRELATION POWER SPECTRUM 

We wish to correlate fluctuations in CO line emission 
with quasar and LRG maps. Our CO temperature mod- 
eling will follow Lll. Lll derived the CO brightness 
temperature by calculating the specific intensity of CO 
emission as a line-of-sight integral of the volume emissiv- 
ity. In Lll, the star formation rate density (SFRD) used, 
given by their Eq. 6, is comparable with the value needed 
to reionize the universe at high z. However, since we are 
interested in lower redshifts than in Lll, the model must 
be modified, as we explain below. 

3.1. Model A: CO Luminosity - Halo Mass 

The basic strategy is to construct a model tha t 
matches three key observational inputs (|Carillil 1201 ID : 
the observed correlation between CO luminosity and far- 
infrared luminosity, the far-infrared luminosity-star for- 
mation rate (SFR) correlation, and the observed SFRD 
of the Universe. In order to predict the spatial fluctua- 
tions in the CO brightness temperature, we need to fur- 
ther connect the star formation rate and host halo mass. 

In comparison to the high redshift z ~ 7 case discussed 
in Lll, this estimate should be on more solid ground in 
several respects. The CO luminosity far-infrared corre- 
lation and the correlation between the far-infrared lumi- 
nosity and star formation rate are measured at z ~ — 3; 
their applicability at higher redshift is questionable. In 
particular, the CO luminosity-SFR correlation may drop 
at high redshift, since the low mass galaxies of interest 
may have low metallicity, as well as an insufficient dust 
abundance to shield CO from dissociating radiation. In 
addition, the increased CMB temperature at high red- 
shift may significantly reduce the contrast between CO 
and the CMB. Furthermore, the overall star formation 
rate density is better determined at z ~ 2 than at z ~ 7. 
On the other hand, the simplistic model adopted in Lll 
to connect star formation rate to halo mass is likely less 
applicable at low redshift, where various feedback effects 
such as photoionization heating, supernovae, and AGN 
feedback have had more time to operate. 

Lets turn to some quantitative estimates. The result 
of combining the CO luminosity far-infrared correlation 
with the far-infrared SFR relation is, in Lll, 



L 



CO(1-0) 



3.2 x 10 4 L P 
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M Q yr-! 
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(2) 



In Lll the authors further assumed a very simple model 
to convert SFR to host halo mass: 

SFR = /*-^— = O.^Mgyr- 1 x 



o.i 



0.17 



10 8 yrs 



M 



(3) 



They took the fraction of baryons converted into stars to 
be /* = 0.1, and a (constant) star formation time scale 
of t s = 10 8 yrs. Star formation was assumed to take 
place with equal efficiency in all halos above some some 
minimum host mass, M s f rjln i n . It was found that the im- 
plied star formation rate density is comparable to the 
critical star formation rate density required to keep the 
Universe ionized, suggesting that this simple presc r iption 
is a reliable estimate. Furthermore, iStark et al.l ()2007l ) 
found that a similar model roughly matches the luminos- 
ity function of Lyman Break Galaxies (LGBs) at z = 6. 
Unfortunately, this simple prescription will not provide 
a good estimate of the SFR at lower redshifts, as we dis- 
cuss. 

We now explain why the estimate for (Tqo) m Lll 
cannot be extrapolated to low redshifts. Taken together, 
Equations[2]and[3]imply that Lco(i-o) oc M 3 / 5 . For sim- 
plicity, we assumed that ico(i-o) instead scaled linearly 
with halo mass after matching the implied CO luminosity 
at M = 1O 8 M . Since massive halos are rare at z ~ 7, 
the CO emissivity in this model is dominated by galax- 
ies in low mass halos close to the minimum CO lumi- 
nous halo mass, M COim i n , and so adopting a linear scaling 
rather than a A/ 3 / 5 power law, has relatively little im- 
pact on the mean CO brightness temperature. At z ~ 7 
the linear scaling leads to only a slight overestimate com- 
pared to the M 3 / 5 power law. Together, these assump- 
tions impl ied L co(1 )(M) = 2.8 x 1O 3 L M/(1O 8 M ). 
Note that iLidz et all (|2011[ ) focused on CO(2-l) assum- 
ing that the J = 2 — 1 and J = 1 — lines have the same 
luminosity, which is very conservative. In the optically 
thick, high temperature limit, the brightness tempera- 
ture would be a factor of 8 larger. At z ~ 7 this more 
than makes up for the possible overestimate from the 
linear scaling. 

Now it is easy to see that applying this blindly at z ~ 2 
may lead to problems because significantly more massive 
halos are no longer rare. As a result, if we assume a 
linear scaling when the true scaling is sub-linear, we will 
significantly overestimate the CO emissivity and bright- 



ness temperature. For example, the above L 



CO(l-O) 



■M 



relation gives L C o(M = lO n M ) = 2.8 x 1O 6 L . How- 
ever if we had used the sublinear scaling, we would have 
obtained L co {M = lO n M ) = 1.8 x lO 5 ^©, which is 
lower by a factor of ~ 15. Although the CO luminosity 
halo mass relation used in Lll seems like a plausible esti- 
mate at high redshifts, where low mass galaxies dominate 
the SFRD and CO emissivity, blindly extrapolating it to 
higher halo masses is problematic. 
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Wc derive a more accurate expression by first adjusting 
the SFR-M prescription to be more suitable at low red- 
shifts, where feedback processes will furt her suppress star 
forma tion in low mass halos. Wc follow IWvithe fc Loebl 
(|2003f) in adopting a halo-mass dependent star forma- 
tion efficiency, as sugges ted by the z ~ observations of 
iKauffmann et a l. (2003). In particular, we assume that 
the fraction of gas turned into stars above M s f rjmin scales 
as M 2 / 3 below some mass Mo, and that the star forma- 
tion efficiency is independent of mass above the charac- 
teristic mass Mq. More specifically we assume 



SFR = /* 



SFR = f* 



Ob M p 
f2 m t s 

n b m 



M 
Mo 



5/3 



-;M > M . 



(4) 



In reality /*, M D and t s likely have some redshift de- 
pendence. Here we fix the characteristic mass scale 
Mo at Mo = 5 x 10 1:L Mf^ (close to the mass scale at 
z ~ 2 in the IWvithe fc Loebl (|2003D model), and nor- 
malize the proportionality constant to match 
the observed SFRD at z ~ 2, p\ = O.1M yr _1 Mpc~ 3 
(jHopkins fc Beacoml [20061) . The mass dependent effi- 
ciency reduces the efficiency of star formation in low mass 
halos, reflecting the impact of feedback processes. It is 
convenient that the SFR scales as M 5 / 3 in this model, 
because this yields a linear scaling of Leo M, although 
with significantly lower normalization than the previous 
relation. (This is strictly true only below Mo but these 
halos dominate the CO emissivity even at low z and so it 
is safe to adopt this scaling for all halo masses here.) This 
means that the only thing that changes in this model is 
the brightness temperature normalization (Too), and not 
the bias and Poisson terms. 

Combining the revised SFR model with Equation [2] wc 
find: 



^co(i-o) — 1 x 10 6 i e 



M 



5 x 10 n M Q 



(5) 



This normalization is a factor of ~ 14 lower than in Lll 
but is likely a better match at z ~ 2 where star formation 
occurs mostly in substantially more massive halos than 
at z ~ 7. 

In order to turn this into an expression for the bright- 
ness temperature we can follow the formulas from Lll, 
but with small modifications. First, we will assume the 
optically thick, high-tcmpcrature limit. This should be 
a good approximation for the J = 1 — and 2 — 1 lines 
at the redshifts of interest. We assume that the duty cy- 
cle of CO luminous activity matches the star formation 
duty cycle, i.e., we assume /duty = t s /t age (z), where we 
fix t s ~ 10 8 yrs, and £ ag e(z) is the age of the Universe. 
The result of this calculation with the normalization of 
Eq.Elis 



(Too) (zj) = 0.65fiK 



/coll (Mso,minj z j) 
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H(z., = 2) 



3.4 x 10 9 yr 

tage( z ) 



1 + Zj 



H(zj) 
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where we have scaled to characteristic values at z ~ 2. 
Note that the redshift scaling is cx (1 + zj) 1 / 2 only in 



ji o.i 




Mco.m,„ = io' m 3 



10'° M_ 



Figure 1. The mean CO brightness temperature in Model A as 
a function of redshift. The solid, dotted, and dashed lines denote 
{Tqo) for M C o,min =10 9 , 10 10 , and 10 11 M o , respectively. Note the 
same mean (Tqq) is assumed for both the CO(1-0) and CO(2-l) 
lines. 



the high redshift limit where H{zj) oc (1 + zj) 3 / 2 . The 
results arc plotted in Figure [TJ In this model the star for- 
mation efficiency declines below the characteristic mass 
scale Mo ~ 5 x 10 11 M Q , which may underestimate the 
higher redshift signal if low mass galaxies form stars effi- 
ciently at high redshift and are CO luminous. However, 
our present data sets provide little constraint at very high 
redshift, and so we are less concerned with the predic- 
tions there. Note that (Tco) in this model is 0.1 — 1/iK 
for all z (for M co>m i n = 10 9 M Q ) and thus will be difficult 
to see in a CMB experiment like WMAP. 

3.2. Model B: CO Luminosity - Star Formation Rate 

We recognize that a weak part of the above estimate is 
the simplistic model connecting SFR and host halo mass. 
In this section we take a more empirical approach to es- 
timating the spatially-averaged CO brightness tempera- 
ture and thereby circumvent the need to connect SFR 
and host halo mass. In particular, here we use recent 
determin ations of the star formation rate (SFR) func- 
tion from lSmit et al.l (|2012f ). These authors use measure- 
ments o f the UV luminos ity function along with the Ken- 
nicutt (jKennicuttl 119981 ) relation, connecting SFR and 
UV luminosity to determine the SFR function, i.e., the 
number density of galaxies with a star formation rate be- 
tween SFR and SFR + dSFR. The observed UV luminos- 
ity is corrected for dust attenuation based on the slope 
of the UV continuum spectra, which gives a luminos- 
ity and redshift-dependent extinction correction. They 
fit Schechtcr functions to the resulting SFR functions 
using new UV luminosity function and extinction mea- 
surements at z ~ 4 — 7, and tabulate results from the 
literature at lower redshift (their Table 3). 
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The Schcchtcr function ( Schcchterl 1 1 9 761) is: 



<S>{SFR)dSFR 



f SFR 
* \SFR* 



cxp 



SFR 



SFR* 



dSFR 



SFR* 



^7) 



and is parameterized by a characteristic SFR, SFR*, a 
characteristic number density, <f>*, and a faint-end slope, 
a. The star formation rate density, p*, can be derived 
by integrating the SFR over the Schechter function. As- 
suming that the Schechter function form is maintained 
to arbitrarily low SFRs then gives: 



p* = c/)*SFR*r[2 + a}. 



(8) 



Here r[2 + a] is a Gamma function. 

We can further combine the SFR Schechter function 
with the observed £<co(i-o) — SFR correlation (Equa- 
tion^, to estimate the co-moving (frequency-integrated) 
CO emissivity. As discussed below, we will assume this 
prescription only for sources above some minimum SFR, 
SFR m i n , and that sources below this critical SFR do not 
contribute appreciably to the CO emissivity. Combin- 
ing Equation [2] and Equation [Jj the resulting co-moving 
emissivity is: 



e co(i-o) — <P*L 



SFR* 
lM Q /yr 



3/5 



1.6, 



SFR 

n 



SFR* 



(9) 



In this equation, Lq denotes the luminosity in the CO(l- 
0) line for a SFR of SFR = lM /yr. We fix this 
at L = 3.2 x 1O 4 L as in Equation [21 and as sug- 
gested by z ~ — 3 CO observations. The factor 
r[a + 1.6, SFR/ SFR*} is an Incomplete Gamma Func- 
tion. This reveals the importance of the sub-linear scal- 
ing of Equation [2] for the CO emissivity: combining the 
sub-linear scaling and the Schechter form for the SFR 
function yields a formally divergent CO emissivity as 
SFR — > for a < —1.6. This may be an artifact of 
extrapolating the Schcchtcr form and/or the sub-linear 
Lco(i-o) scaling with SFR to arbitrarily low SFR; pre- 
sumably one or both of these relations drop-off at low 
SFR. For instance, low luminosity dwarf galaxies may 
have small metallicites and fall below the extrapolation 
of the Leo — SFR relation . These relations generalize the 
approach of iCarillil (|2011| ). who estimated the CO emis- 
sivity from the SFRD assuming the two scale linearly 
with SFR. Our more accurate relation requires, however, 
specifying SFR m i n to relate the co-moving CO emissivity 
and the SFRD. 

Using the Equations in Lll we can relate the co- 
moving CO emissivity as calculated above to the 
spatially-averaged CO brightness temperature. It is use- 
ful to note that this gives: 



(Too) (zj) = 0.65/iA' 



e CO(l-0) 



H(zj 



6.3 x lO^oMpc -3 

2 



H(zj) 



1 + zj 



(10) 



We calculate the results of Equations [9] and [TU] from 
z = 0.2 — 6.8 us ing the best-fit results in Table 3 of 
iSmit et all (|20l2l ). We focus on the UV- and MIR- 
derived SFR functions from that paper, and do not in- 
clude the Ha-based estimates, although these appear 
to be comparable. We ignore the point at z = in 
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SFR„ 

/ SFR m 

/ SFR m 

SFR m 



= 0.01 M s /yr 
= 0.1 M /yr 
= 1.0 M e /yr 
= 5.0 M /yr 



Figure 2. The average CO brightness temper atur e in Model B. 
These results are estimated fr o m Eq uations[9land ll0l using the SFR 
function fits from Smit et al. (20f2). The solid pentagon at z = 
is the average CO brightness temperatur e inferred from th e CO 
luminosity function measured at z = by[Kcres et al. 1 120031) . 



their Table 3 from iBothwell et al.l (|2011l) , since there ap- 
pears to be a typo in the fit parameters from this pa- 
per, as remarked in the footnote to Table 3. The results 
of these calculations are shown in Fig. [2] The overall 
rcdshift evolution shown here is noisy because (Tco) is 
estimated from several discrete observational measure- 
ments, each with significant observational error bars. 
As anticipated above, the brightness temperature esti- 
mates arc sensitive to SFRmin, which is somewhat uncer- 
tain. The observed Lco(i-o) — SFR correlation includes 
galax ies with SFRs down to ~ O.5M /yr (jWang et al.1 
[2O1O0, but the data are a bit sparse at low SFR. The 
Schechter function fits probe galaxies down to typical 
SFRs of ~ 0.1 — lM /yr, with some reds hift depen- 
dence at the faint-end limit (e.g. Fig. 2 of ISmit et al.l 
(|2012t) ) . The results are particularly sensitive to SFRmin 
aX z = 6.8, where the slope of the faint-end luminosity 
function is especially steep (a = —1.96). The more em- 
pirical model shown here agrees with the models of the 
previous section of Figure [TJ at the order-of-magnitude 
level, although the results in the SFR m i n = O.OlM /yr 
case are larger by a factor of ~ 4 at z ~ 2 than the 
earlier Mcomin = lO 9 A/ /yr model. The high-redshift 
results at z > 5 are also broadly consistent with the es- 
timates in Lll. There are still significant uncertainties 
in the Schechter function fits to the SFRs, and in the 
extrapolations to the faint end, but this latter estimate 
is probably more secure than the estimates of the pre- 
vious section which rely on an oversimplified model to 
connect SFR and halo mass. Note that we will estimate 
power spectra for Model B by multiplying the Model 

A spectra by (T co ) 2 (Model B)/ (T co ) 2 (Model A), al- 
though strictly speaking this is not accurate since the 
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clustering bias depends on the SFR-M relation, which is 
not determined for Model B. 

As one final sanity check, we can us e the mea- 
sured CO(l-O) luminosity function from iKeres et al.1 
at z — to estimate the CO emissivity and 
brightness temperature at z = 0. They provide a 
Schechter function fit with p(L) = d$/d\og 10 L = 
p*MW)(L/L*) a+1 cxp(-L/ L*), and = 7.2 x lO" 4 
Mpc~ 3 mag" 1 , a = -1.30. IKeres et all (|2003f ) quotes 
results in terms of velocity-integrated CO luminosities, 
with a characteristic value of Lco.v-.* — 1.0 x 10 7 Jy 
km s _1 Mpc 2 . We convert thi s into solar units u s ing th e 
relations in the Appendix of lObrcschko w et al.l (|2009f ). 
finding = 9.6 x 10 4 L f 7 1 . Using the best-fi t Schechter 
function parameters from IKeres et al . (2003), and inte- 
grating to £co(i-o) = gives (Tco) = 0.064p,K, broadly 
consistent with our estimates, extrapolated to z = (see 
the solid pentagon in Fig. [2]). 

3.3. CO Clustering 

In order to model the clustering of CO emitters, fol- 
lowing Lll, we will assume the standard relation between 
CO luminosity and host halo mass, linear bias, that the 
scales of interest are much larger than the virial radius of 
the relevant halos so that we can neglect the non-linear 
term, and that the halo shot noise obeys Poisson statis- 
tics. With this assumption, the 3D power spectra for CO 
temperature fluctuations goes as 



Pco(k,z) = (T C of(z) 



b 2 {z)P lin (k,z) + 



1 



(M 2 ) 



/duty {z) (M) 2 
(11 



where b(z) is the effective z-dependent halo bias given in 
Eq. 15 of Lll 



b(z) = 



IZ,^ dMM inKM,z) 



M,., 



dMM^ 



(12) 



where dn/dM is the halo mass function from 
iTinker et al.l (|2008|) and M is the associated halo mass, 
and Pij n (fc, z) is the linear matter power spectrum. The 
second term in brackets for Pco is the shot noise, with 



<M 2 > = 



(M)- 



dM NT 



dM M 



dn 
dM ' 

dn 
dM ' 



(13) 



We implicitly assume that every dark matter halo will 
host at least one CO emitter with a duty cycle /duty(z). 
We plot the CO three-dimensional power spectrum in 
Fig. |31 We will also assume linear biasing for quasars so 
that the 3D quasar and LRG spectra are 



P Q (k,z) 

PLRG(k,z) 



■ b 2 Q {z)F\ ia (k,z) 

' ^IflG^lin^ z ) ? 



(14) 



where &q(z) is the quasar clusterin g bias that w e com - 
pute using the estimate of Eq. 15 in lCroom et al.l ([2005D 
and bhRG is the LRG clustering bias. In general, the 
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Figure 3. The 3D power spectrum for CO fluctuations in the 
rcdshift range < z < 4. The top panel is for Model A with 
-Mco.min = 10 9 Mq, while the bottom panel is for Model B with 
SFR mln = 0.01M o /yr. The black dashed (dotted) line shows the 
clustering (shot noise) term for 2 = 2. 



LRG bias would be redshift-dependent, but we will as- 
sume a constant bias for the LRGs. Similarly, the 3D 
cross power-spectra can be written as 

Pco-Q(k, z) = r Q (T C o) (z)b(z)b Q (z)P Vm (k, z) 
Pco~LRG(k,z) = r LRG (Tqo) (z)b(z)b LRG P Un (k, z\\§) 

where we include cross-correlation coefficients tq and 
tlrg between CO emitters and quasars and LRGs, re- 
spectively, which we assume to be scale- and redshift- 
independent throughout this analysis. Note it is possi- 
ble that quasars and LRGs may live in the same halos 
as CO emitters, contributing a shot noise term to their 
cross-correlation. We will neglect this shot noise in the 
analysis. 

Assuming the Limber approximation in the small scale 
limit (typically £ > 10), we can then derive the angular 
auto/cross power spectra. The cross-power spectra have 
the form 



CO- 



= 1 dz mM$m x 



c x 2 (z) 
P co _ [k = £/ X (z),z\ 



C 



CO-LRG _ 



= dz 



H(z) fco(z)f L RG(z) 
c x 2 (z) 
Pco-LRc[k = £/x(z), z] , 



(16) 
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where H(z) is the Hubble parameter and /co(z), Iq{ z )i 
and Jlrg(z) are selection functions of the CO tempera- 
ture fluctuations, the quasar distribution, and the LRG 
distribution, respectively. In our analysis, we assume 
flat, normalized selection functions for CO, quasars, and 
LRGs as a rough estimate. The three relevant angular 
auto-power spectra are given by 



nCO _ 
— 

qLRG _ 



dz Hiz)f ¥ z) Pco[k = £/ X (z),z] 



dz 



c x 2 0) 
H{z) fU z ) 
c x 2 (z) 



P Q [k=£/ X (z),< 



dz H{z) P LRG [k = t/ X (z),z] (17) 



We plot the predicted angular power cross-spectra for 
Model A with quasars in Fig.|4]and the cross-spectra with 
LRGs for Mco.min = 10 9 Af o in Fig. [5] On the plot we 
superimposed an estimate for Model B by just rescaling 
(Tco)- For the LRGs we assume b^ua = 2.48, which we 
measured from the data (see Sec. I5.2|) . Given the finite 
redshift widths defined by the WMAP bands, the evolu- 
tion of the peaks of the angular spectra with the band 
frequency reflects the evolution in angular diameter dis- 
tance with redshift. A higher band frequency means a 
lower redshift, which means a larger angular scale (lower 
£) for the peak. For each band, the CO(l-O) spectrum 
can be higher or lower than the CO (2-1) spectrum de- 
pending on the redshift ranges probed, which determine 
the halo, quasar, and LRG biases and the CO brightness 
temperature. 

Given Model A, we also find that 



C 



CO(l-O) 



(r CO{2- 
\°£=200 



£=200 

), 3.48xlCT 7 
and 7.80 xlCT 8 



3.90xl(T 7 (2.64x10 
(3.88xl0" 7 ), 1.93xl0- 7 (3.73xl0~ 7 ) 
(2.83xl0^ 7 ) ^K 2 respectively in the Ka, Q, V and W 
bands if we neglect the shot noise term, which becomes 
dominant only at £ = 579 (N/A), 292 (N/A), all £ (744), 
and all £ (170) for the same bands. Note that these 
numbers change if M COjm i n ^ 10 9 M Q . For reference, 
the CMB angular power spectra is Cf MB ~ 1 /iK 2 at 
£=200 so it will act as an important source of noise, 
if not subtracted. High quasar shot noise causes CO 
and quasars to not be perfectly correlated even on the 
largest scales. Wc note that the shot noise is more 
important at lower z (lower v) as the number of massive 
halos increase more rapidly than the number of small 
mass halos. As will be discussed below, the instrumental 
noise at £ = 200 is 0.0253, 0.0230, 0.0347, and 0.0501 
/iK 2 . In Figs. E] and El we also present the various 
sources of noise in comparison to the cross-correlation 
signal. Note that we do not include foreground noise, 
although it is included in the final results. Other than 
foreground noise, we see that the spectra are dominated 
by the WMAP instrumental noise. 

Given these numbers, the prospect of measuring a 
cross-correlation signal are not very high given the cur- 
rent state of our model. However, given the theoretical 
uncertainties, we will still attempt to measure directly 
this correlation with currently available data. 

4. MEASUREMENT METHODOLOGY 
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Figure 4. The predicted angular cross-power spectra and noise 
curves with each WMAP band cross-correlated with its correspond- 
ing quasar map. The solid, red line is the cross-correlation sig- 
nal for CO line emission for Model A with M co m i n = 10 9 Mq, 
and the purple, short-dashed line is for Model B with SFRmin = 
O.OlA#0/yr. The blue lines are the noise curves assuming Model A 
for the following cases: cosmic variance and CO shot noise limited 
(dotted), including WMAP instrumental noise (dot-dashed) only, 
including quasar shot noise only (dot-dot-dot-dashed), and includ- 
ing both noise sources (dashed). Note that the noise curves are 
slightly modified for Model B. 



We estimate the CO-LSS tracer angular power spec- 
trum, where the LSS tracer can be either quasars or 
LRGs, using a minimum- variance estimator of the form 



C 



CO-Tr 



E 

m——l 



n co n Tr * 



(2£ + l)f sky , 



CO-Tr 



(18) 



where J s k y ,co-Q = 0.140 and / s k y ,co-LRG = 0.197 are 
the sky coverage fraction of the WMAPxSDSS map inter- 
sections for quasars and LRGs, respectively, and af m are 
the invcrsc-noisc-filtcrcd spherical harmonic coefficients 
for observable X , given by 



-x 



(Cf +C^ x )[(S x + N x r 1 d x ] 



X\lm , 



(19) 



where and C™' X are the fiducial signal and noise 
angular auto-power spectra, respectively. We will as- 
sume Model A in Sec. O with M COlIoin = 10 9 M Q for 
the fiducial signal auto-power spectrum, noting that the 
estimator should not depend significantly on the fiducial 
spectrum. The CO thermal noise and quasar shot noise 
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Figure 5. The predicted angular cross-power spectrum and noise 
curves for Model A with the WMAP W band cross-correlated with 
the LRG map for M co m ; n = 10 9 Mq. The solid, red line is the 
cross-correlation signal for CO line emission, and the purple, short- 
dashed line is for Model B with SFR min = 0.01Af Q /yr. The blue 
lines are the noise curves for the following cases: cosmic variance 
and CO shot noise limited (dotted), including WMAP instrumental 
noise (dot-dashed) only, including quasar shot noise only (dot-dot- 
dot-dashed), and including all noise (dashed). Note that the noise 
curves would be slightly modified for Model B. 



angular power spectra are 



C 



n,COi 



and 



C 



n.Trt 



(20) 



(21) 



for the i band, where j is the band being used as the 
CMB map (see subsequent paragraph). Af2 is the pixel 
size, a\. . is the average of a\, + a\. over all pixels, n, tr 
is the average number of objects per steradian for the 
tracer, and Wj l includes the pixel and WMAP beam 
window functions, dx is the data vector for observable 
X with the entries in dco and d^r being the CO fluctu- 
ation (in fiK units) and the object number overdensity 
for the tracer, respectively. Sx is the fiducial signal co- 
variance matrix and Nj is the noise covariancc matrix 
with (d x d£) = S x + N x - (Sx + N x ) _1 works as an 
operator that weights each data point by its covariance, 
effectively filtering out very noisy modes. It is evalu- 
at ed using a multigri d preconditioner (see Appendix A 
of ISmith et all (120071) for details). The variance for the 
estimator in Eq. [T5] is given by 



Var^ 



where 



CO-Qi 



1 



[2i + l)/ sky 



CO-Q-,2 



r + 



D?°D? 



22) 



D 



CO 



E 



= CO|2 



Dl r = E 



^(2^+l)/ sky , co 

I^Tr 12 



L im\ 



^ ( 2 £+l)f skytTl 

ra— — l 



(23) 



and /sky, co = 0.783, / s k y ,Q = 0.17, and / s k y ,LRG = 0.203 
are the sky coverage fractions of our WMAP and SDSS 
data samples, respectively. For completeness, we also 



Var[£>f °] : 
Vax[bJ r ] 



(2£+l)/ skv , co 
2 

(2^+l)/ sky ,T^ 



(24) 



By estimating the errors this way, we include all sources 
of noise including galactic foregrounds. 

As discussed above, the WMAP temperature maps arc 
dominated by the CMB and galactic foregrounds, with 
CO emission constituting a small contribution. CMB 
fluctuations would drastically increase statistical errors 
in our CO search, so we choose to try and subtract out 
the CMB. Because the WMAP bands are so large, we 
cannot model the CMB fluctuations to subtract them 
properly. However, since the V and W maps are dom- 
inated by the CMB outside the galactic masl<0, we can 
use these maps to subtract the CMB from all the maps. 
Specifically, we subtract the W map from the maps Ka 
and V, and we subtract the V map from the Q and W 
maps. Note that all the maps have zero mean. Since 
the WMAP maps are already given in physical tempera- 
tures, we can just subtract them directly before convert- 
ing them to brightness temperatures (see Sec. I2.1[) . The 
noise fluctuations in the V and W maps make the sub- 
tractions imperfect, even if we were to assume the fore- 
grounds in these maps were negligible. This causes noise 
fluctuations to increase due to the noise introduced by 
the subtracted map according to Eq. but this should 
not introduce a bias. However, subtracting different 
bands, even with the extra noise, is worthwhile because 
keeping the CMB perturbations would contribute much 
more noise. We give an example of a CMB-subtracted 
map in Fig. [5] From the two maps, it is evident that 
our method removes the hot and cold spots of the CMB. 
Note that this subtraction does not affect the CO signal 
we aim at cross-correlating with since the bands do not 
overlap in CO redshift sensitivity as is clear in Tab. [TJ 

5. RESULTS 

5.1. Quasar Results 

We present measurements of the CO-Q cross- 
correlation in Fig. [7] Remember that the Ka band mea- 
surement for the CO (2-1) line was not performed. Even 
though the CMB was subtracted from the data, the error 
bars on the cross-correlation are still much larger than 
the predicted spectra from Model A. The WMAP tem- 
perature noise and resolution, as well as the number den- 
sity of quasars and any remaining fluctuations due to 
foregrounds, affect the statistical errors, but from Fig. @] 
we can infer the noise from the brightness temperature 
residuals is by far the greatest contributor compared with 
its theoretical noise-free fluctuations. It appears that our 
measured angular power spectra are consistent with a 
null signal and Model A in all the redshift bins. Model B 
is fairly high in some of the bands, but not entirely ruled 
out. We also show lOx the Model B estimate, which we 
see is ruled out in several of the bands on the larger scales. 
This becomes much clearer in the (Too) constraints be- 
low. Thus we can infer that the CO brightness is not 

4 Specifically, the foreground contribution to the V and W maps 
are approximately 20%. 
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Figure 6. The brightness temperature map in the WMAP Q band 
(left) and the same map with the brightness temperature map from 
the W band (estimating the CMB fluctuations) subtracted (right) 
with galactic emission and point sources masked out. The color 
legend is given in units of mK. 

significantly brighter than Model B (e.g. there is no abun- 
dant population of faint CO emitters that is not included 
in our model). However, a future experiment with an 
increased temperature sensitivity and quasar density is 
needed to detect the CO-Q cross-correlation (see Sec. [5]). 

Using the extra simplifying assumption that (Too) is 
constant over each redshift window we can translate our 
measurement into a constraint on the effective (Woo) 
temperature, (bTco)- We define the (bTco) estimator 
by rewriting Eq. [16] as 



n CO-Q 



'(bTco) 



dzm f C o( Z )f Q ( Z ) bQ{z) x 



X 2 (z) 



Pnn(k = £/ X (z),z) 



= r(bT co )F e . (25) 

Assuming r = 1, we then can construct a minimum- 
variance estimator for (bTco) according to the expres- 
sion 



(bTco) = 
with uncertainty 



£,(7f°- Q iVVar[Cf°- 



1 



'bT co 



£ £ F//Var[G- 

= E 



co-c 



Ft 

1 u 



Var[C; 



CO- 



(26) 



(27) 



In practice, we limit ourselves to 80 < £ < 1000 to avoid 
the Integrate d Sachs- Wolfe (ISW ) "contamination" on 
larger scales (|Sachs fc W olfe 1963). Using the same for- 
malism as above, our constraints on (bTco) (assuming 
r = 1) are given in Table [2j If we further assume knowl- 
edge of the bias, which in practice means that we know 
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Figure 7. The measured, binned CO-Q cross-correlation angular 
power spectrum measurements with 1-sigma error bars. The solid 
line is the nonzero theoretical CO-Q correlation from Model A 
for M co m i n = 10 9 Mq, the dotted line is the rough estimate for 



Model B for 5Fi? min = 
estimate for Model B. 



0.01MQ/yr, and the dashed lino is lOx the 



that host halo mass of CO emitting objects, these con- 
straints yield Too constraints, where we rewrite Eq. [T5] 

as 



rtCO-Q 7fi 

C v = rT C o 



, H(z) fco(z)f Q (z) 

az 5—- b{z)bQ(z)x 

c x (z) 



Pi in (k = e/ X (z),z). (28) 

Still assuming r = 1, we plot constraints of (Too) 
for each line over four redshift bins, comparing the re- 
sult to Model A for M CO)inin = 1O 9 M , Model B for 
SFR min = 0.01A:T Q /yr, in Fig. H and lOx Model B. It 
is evident from these plots that the (Too) signal from 
Model A and Model B are not detectable from our anal- 
ysis, but we know that the brightness cannot be an or- 
der of magnitude greater than Model B. In fact, a model 
more than 3 times greater than Model B would be ruled 
out by the combined constraints of the two higher red- 
shift points for the CO(l-O) line. Another way of saying 
this is that the signal-to-noise ratio (SNR) for the two 
points, which is the SNRs of the two individual points 
added in quadrature, where the signal is the difference 
between the data point and the model, is more than 3 
for a model more than 3 times Model B. 

5.2. LRG Results 
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Table 2 

Measured bTco with 1-sigma error bars of CO emission lines 
J = 1 — > and J = 2 — > 1 from cross-correlations of WMAP 
temperature bands with SDSS DR6 quasar (QSO) counts. 



WMAP band 


bTcoa-O) (mK) 


bT C 0(2-i) (MK) 


Ka 


3.89 ± 19.39 


N/A 


Q 


-11.5 ±9.8 


-138. ± 207. 


V 


1.41 ±8.68 


18.7 ±24.1 


w 


2.98 ±3.94 


-0.74 ±6.36 
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Figure 8. Limits on (T C q) of CO(1-0) and CO(2-l) over four 
redshift bins. The solid line is the fiducial temperature from Model 
A for M co . 

min — lO^il^f©, the dotted line is the fiducial temperature 



O.OlAf0/yr, and the dashed line is 



from Model B for SFR n 
lOx the estimate for Model B. 



Before measuring CO parameters using the LRG cor- 
relation, we proceed to determine the LRG clustering 
bias bLRG for the full sample, ignoring the bias redshift 
evolution. We start by rewriting the LRG angular auto- 
spectrum in Eq. [T71 as 



r<LRG _ i2 

°£ ~ U LRG 
— U LRG 



d 

riLRG 



H(z) fl RG (z 
c x 2 0) 



-Pun[k = e/ X (z),z] 



(29) 



where we assume a constant LRG bias. We use a 
model for the measured LRG auto-spectrum given by 

Df RG = b 2 LRG Gf RG + tf RGshot , where tf RGshot is 
the shot noise component for LRGs given in Eq. [21] This 
model allows us to construct the estimator for the LRG 
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Figure 9. The measured, binned CO-LRG cross-correlation an- 
gular power spectrum measurements with error bars. The solid 
line is the nonzero theoretical CO-Q correlation from Model A for 
A'/ C o,mm = 10 9 A'/q, the dotted lino is the rough estimate for Model 
B for SFii m i n = O.OlM0/yr, and the dashed line is lOx the esti- 
mate for Model B. 



bias 



[b ]LRG = sw^n — ^ p) 



with uncertainty 



12 51 



{G LR G) 2 



Wlbg t Var[^« G ] ' 



(31) 



where we limit the sum to 10 < £ < 100. Using this 
method, we find bLRG = 2.4820.055, which is a little 
high compared to other analyses bu t still cons i stent given 
that our redshift range is wider. In lHo et al.l (|2008| ). the 
measured LRG bias for the relevant redshift range was 
bLRG = 2.03±0.07. We list results assuming both values. 

With an LRG bias, we can perform the CO-LRG cross- 
correlation measurement, the results of which we show 
in Fig. [5J These measurements, like the CO-quasar mea- 
surements, are consistent with both models, as well as 
lOx Model B, making it not really useful. We estimate 
CO parameters for 0.16 < z < 0.47 in a similar man- 
ner as the previous section. The values we measured 
for bT CO ( 1 _Q\ and T C q^i_^ assuming our LRG bias are 
0.76±1.06 /^K and 0.56±0.78 llK, r espectively. Assuming 
the LRG bias from lHo et al.l ( 20081 ) . the values change to 
0.93 ± 1.30 fiK and 0.69 ± 0.96 [iK, respectively. Note 
that the cross-correlation measurement was performed 
without including the measured LRG bias to the fiducial 
signal; however, we repeated the measurement with the 
bias and confirmed that the final result was unchanged. 
These constraints are indeed tighter than those from the 
quasars, due to a much higher areal density of objects. 
However, at z ~ 0.25, the relevant redshift, the predicted 
(Too) for both models are in the 0.1-0.2 range, making 
these constraints still not useful. 

5.3. Possible Contaminants 

Our measurements could be contaminated by any ex- 
tragalactic foregrounds correlated with quasars or LRGs. 
We already attempted to remove ISW by neglecting 
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large-scale modes, but other foregrounds may contam- 
inate our signal. For example, quasars could live in SZ- 
contributing clusters or be closely correlated with dust- 
rich galaxies, which would lead to a bias in our mea- 
surements. Since we subtract maps, it is possible for a 
negative bias to occur, canceling a true signal. Although 
this is unlikely, if nothing else it is useful to describe these 
possible contaminants since they will be important in fu- 
ture analyses. We will briefly discuss these possibilities 
below. 

As mentioned in S ec. 12. 1| the thermal Su nyaev- 
Zel'dovich effect (tSZ) (jSunvaev fc ZeldovichllT980l ) is ex- 
pected to correlate CMB temperature with large scale 
structure tracers like quasars or LRGs. In fact, there are 
hints of a SZ-quasar cross-c orrelation in WMAP/SDSS 
data (fChattcrj ee et al.ll2010P) . Thus, it is important to 
check the level of tSZ contamination in our measurement. 
Since tSZ has a well-defined frequency dependence, we 
can rule out tSZ if our cross-correlation measurements do 
not follow this frequency dependence. Furthermore, we 
note that while performing the CMB subtraction using 
the V or W band, we also alter the tSZ signal frequency 
dependence so that at every frequency, X, we will write 
the tSZ amplitude, as AT SZ = ATg| - AT^ or w . Even 
though we do not expect tSZ in our CO analysis to be 
significant, in the following analysis we confirm whether 
or not this effect is indeed a problem. 

The tSZ effect causes the CMB temperature to receive 
a secondary perturbation ATsz every time it scatters 
with an object, with the perturbation given by 



AT sz (x) = yT CMB f{x) , 



(32) 



where y is the object-dependent Compton y parameter, 
x = hv/kT, and f(x) — xcothx — 4. This implies that 
once the CMB reaches us, its tSZ perturbation in a pixel 
due to a set of objects (of one type) in a redshift bin will 
be 5T = ATsz (A — N), where A is the number of objects 
in the pixel and ./V is the average number of objects per 
pixel. This model allows us to relate the angular cross- 
correlation between tSZ and a LSS tracer, Tr, to the 
tracer's auto-spectrum according to the form 



C: 



SZ-Tr 



Af sz N(C, 



Tr 



£jTr, shot ^ 



(33) 



Using this model, we can construct an estimator for ATgz 
using quasars as the LSS tracer, given by 



AT< 



sz 



SZ-Q /r ,Q 



{C? + C £ Q ' shot )/Var[C J 



A 



Y.t(Pt + Cf shot ) 2 /Var[C^- y ] 



SZ-Q-, 



(34) 



The procedure for measuring the tSZ-Q cross-spectrum 
is equivalent to measuring the CO-Q cross-spectrum, ex- 
cept that tSZ-Q uses the CMB physical temperature 
while CO-Q uses the brightness temperature. This al- 

lows us to set Cp = C f /fbi, where /b r is the 
prefactor in Eq. [TJ for converting physical temperatures 
to brightness temperatures. 

We show the measured ATsz at each redshift for the 
quasar samples in Fig. QJJJ For the LRG sample, we find 
ATsz = {0.54±1.05,0.61±1.12}for b LRG = {2.48,2.03}, 
respectively. None of these measurements are significant, 
but since we did not detect CO emission, we can't rule 
out the presence of tSZ either. We emphasize that it is 



CO(1-0) Sample 




C0(2-l) Sample 




Figure 10. Limits on AT SZ for the CO(l-O) and CO(2-l) quasar 
samples. For the CO(l-O) samples, the points from left to right 
to right are for the WMAP maps W-V, V-W, Q-W, and Ka-W. 
For the CO(2-l) samples, the points from left to right to right arc 
for the WMAP maps W-V, V-W, and Q-W. The z axis labels the 
redshift of the quasar map. 

necessary to search for tSZ if CO line emission is detected 
in any future experiment. We note that a spectrograph 
like the one we propose in Sec. [S] will have a high enough 
resolution to easily remove tSZ (and the CMB) directly, 
so this should not be an issue. 

We also investigated contamination by dusty galaxies. 
Dust has a positive spectral index, causing its contri- 
bution with the W and V bands to be higher than the 
other bands. Thus, while correlating the difference of 
two WMAP bands, e.g. Q-V, with quasars, we may get 
a negative correlation which could possibly cancel the 
CO-quasar correlation. We test this possibility by cor- 
relating Ka-W, Q-V, and V-W with quasars that would 
correlate with the CO(l-O) line in the Q, Ka, and Q 
bands, respectively. Since this correlation cannot come 
from CO, we know that if dusty galaxies are canceling 
the CO-quasar correlation, they would create a negative 
signal. Note that W-V correlated with quasars or LRGs 
should not have this problem because the signal would 
be positive. We find the correlations are indeed consis- 
tent with zero, as seen in Fig. [Til so we can conclude 
that dusty galaxies are not canceling a CO-quasar cross- 
correlation detection. Similarly, this null result also rules 
out other possible contaminants such as radio sources. 

5.4. Another use of WMAP data 
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Figure 11. The measured, binned CO-Q cross-correlation an- 
gular power spectrum measurements for out-of-band quasars with 
1-sigma error bars. This measurement shows that the tempera- 
ture correlation with dust is not significant enough to bias the 
CO-emission measurements. The dashed line represents a null cor- 
relation. 



Another method of detecting CO temperature fluctu- 
ations is by cross-correlating the CO(l-O) and CO(2-l) 
lines coming from the same emission redshift. For each 
WMAP band, the CO(l-O) redshift range and the CO(2- 
1) redshift range arc different; this can cause the CO(l-O) 
redshift range for one band to overlap with the CO (2- 
1) redshift range of another band, e.g. Ka for CO(l- 
0) and V for CO(2-l). Since these two rotational lines 
would presumably be emitted by the same sources, we 
can cross-correlate pairs of WMAP bands that overlap in 
this way to search for CO emission. This could be a way 
to measure the line ratio as a function of redshift. Specif- 
ically, we could search for an angular cross-correlation of 
the form 



C\ 



CO-CO' 



, H{z) fco(z)fco'(z) 
dz — x 

c x (z) 

sPco[k = £/x(z),z], 



(35) 

integrated over the intersection of the two redshift ranges 
where unprimed CO is CO(l-O) and primed CO is CO(2- 
1). We include a parameter s, the ratio between the 
two CO lines. We can attempt to measure a unitless 

amplitude B such that Cf ' 00 ' = BC^~ co ' . In this 
case, B has an estimator similar to Eqs.l2"6landl2"Tl In this 
procedure, there would be other sources of correlation in 
the WMAP bands coming from various foregrounds, with 
the CMB as the dominant foreground. To remove this 
complication, we can either remove these foregrounds by 
hand or cross-correlate pairs of WMAP bands that do 
not correspond to the associated CO redshift window as 
a cross-check to see if the difference in cross-correlation 
between the two cases is statistically significant. A third 
option would be to use a three point function, something 
like Ka-V-QSO, as this would lead to an extra handle on 
the foreground correlation. 

Based on the WMAP band CO redshift ranges in Ta- 
blc[TJ we find that the Ka band for CO(l-O) should over- 
lap with the V band for CO(2-l), while the Q band for 
CO(l-O) should overlap with the W band for CO(2-l). 
The other combinations (Ka-Q, Ka-W, Q-V, and V-W) 
should not exhibit any cross-correlations that are due to 
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Figure 12. The predicted angular cross-power spectra from 
Model A for M COjmin = 10 9 M o with a WMAP band for CO(l- 
0) cross-correlated with another WMAP band for CO(2-l). The 
solid line is the clustering term and the dotted line is the shot noise 
term. 

common CO emission. The predicted spectra for Ka-V 
and Q-W are plotted in Fig. [T3J We attempt to subtract 
the CMB from all four bands, cross-correlate the Ka-V 
and Q-W pairs, and then look for a detection of an am- 
plitude. However, when we calculate the variance of the 
cross-correlation Cf °~ co given by 



CO-CO' 



{21 + l)/ sky ,co 



(C 



CO — CO' \2 



D?°D?°' , (36) 



we find extremely large errors for this CO amplitude. 
The errors for the amplitude for Model A are ub — 1400 
for Ka-V and erg ~ 1200 for Q-W. The numbers are 
much better for Model B, being <jb — 45 for Ka-V and 
ob — 68 for Q-W, although still less powerful than those 
from the WMAP-quasar cross-correlation. The limits 
from the COxQSO analysis greatly outperform the lim- 
its from a potential WMAP CO(l-0)xCO(2-l) analysis 
because the instrumental noise is relatively much higher 
than the shot noise as shown in Fig. |4j In addition, 
WMAP foregrounds and CO emission from higher J can 
contaminate the signal, so we do not attempt to measure 
it in this analysis. 

6. CO EXPERIMENT FORECASTS 

Now that we have determined WMAP's limited ability 
to constrain CO temperature, we investigate what can 
be done in future experiments. We know that the lim- 
iting factor in the CO-Q cross-correlation measurement 
is how well we can reduce temperature noise as well as 
increase the quasar density. Unlike WMAP, a dedicated 
experiment should have small frequency bands in order 
to maximize power from the CO line and to subtract 
continuum sources, including the CMB and other fore- 
grounds. We consider such a configuration here- 
in this section, we seek to determine the signal-to-noise 
ratio (SNR) of the cross-correlation between the bright- 
ness temperature fluctuation of CO(TO) measured by a 
hypothetical CO experiment and the latest spectroscopic 

5 Mike Sciffcrt, private communication 
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Table 


3 




CO experiment 


parameters. 


Survey Parameter 


Expt 1 


Expt 2 


CT,fwhm 


6.7 


4.7 


Sfwhm (arcmin) 


40.5 


28.3 


Survey Area (deg 2 ) 


550 


270 



LSS surveys at z ~ 3. This CO experiment measures 
the CO(l-O) rotational line at observed frequency 28.8 
GHz with a 1 GHz bandwidth and 20 50-MHz frequency 
channels. At a frequency resolution of approximately 
R = 600, this experiment should be able to effectively 
model CMB and foreground emission and remove their 
contributions. The central frequency of the band would 
correspond to z = 3 with channel widths corresponding 
to Az = 0.00694 (2.9306 < z < 3.0694). Table shows 
sample parameters of this experiment for two cases. For 
each case, we set a instrumental error CT.fwhm, a beam 
FWHM, and a survey area. Note that 0T,fwhm i s f° r a 
pixel corresponding to the FWHM we choose. We then 
have a noise angular power spectra of the form 

C7' co = 4 fwhm (O.42450 fwhm ) 2 /Hf . (37) 

We begin with the BOSS final spectroscopic quasar 
survey. The BOSS full survey will cover 10,200 deg 2 
by 2014, long before the time when the CO experiment 
would start observing. In the range compatible with our 
hypothetical experiment, based on projections from the 
BOSS Year One (partial) survey, we expect BOSS to 
detect 5465 spectro-Q SOs per steradian . Based on the 
redshift distribution in lRoss et al.l (|2012l ). we see that the 
distribution is pretty flat in our range of interest, making 
us set the density in each redshift bin equal to 5465/20 
= 273 str" 1 . 

In principle, the spectral resolution is such that the 
cross-correlation will be performed directly in three di- 
mensions, as is currently performed w ith 21cm surveys 
(jChang et alj|2010t iMasui et alj|2012j ) and Lll. For the 
sake of simplicity, we will however compute here a sim- 
ple forecast using a bin-by-bin 2D angular power spectra. 
Using the 2D angular power spectrum to constrain CO 
temperature neglects modes along the line of sight, low- 
ering the SNR. Thus, we take the SNR forecasts in this 
section to be lower limits to what could be achieved with 
a full 3D power spectrum analysis. 

The SNR for the experiment is then 

(c co-Q^ 2 

SNR2 = ^-binsE 4r[cr°-V (38) 

where iV z _bins = 20 and we assume each bin contributes 
approximately equally. Note that we neglect the cross- 
correlation shot noise in this calculation, as well as for the 
rest of this Section. We find for Model A with M co , min = 
10 9 M Q the forecasts SNR = 1.2 (5.4) for Experiment 1 
and SNR =1.8 (8.0) for Experiment 2 for each redshift 
bin (over the full redshift range). Experiment 2 may 
be of interest for the full redshift range, but systematics 
resulting from not subtracting SZ and dust contamina- 
tion properly may degrade the signal. An autocorrelation 
would have an even lower SNR (SNR ~ 7.1 for Experi- 
ment 2). Probing Model B should have better prospects; 



although we can't estimate the cross-spectrum directly, 
we can estimate that since (Tco) at z = 3 for Model 
B with SFR min = 0.01M©/yr is about 4/0.85 ~ 4.8 
higher than for Model A with M COimin = 1O 9 M , its 
cross-spectrum will be approximately 4.8 times higher. 
This changes the SNR values to 4.2 (18.8) for Experi- 
ment 1 and 4.4 (19.7) for Experiment 2. Note that the 
SNR does not increase by a factor of 4.8 because the noise 
is model-dependent. The autocorrelation SNR would ac- 
tually be higher than that for the cross-spectra in this 
case (SNR ~ 49 for Experiment 2), but foregrounds can 
degrade this signal. Another option would be to have 
an instrument that probed CO fluctuations at two fre- 
quencies with one twice the frequency of the other. With 
this, you could search for a cross-correlation signal from 
CO(l-0)xCO(2-l); the SNR for an "equivalent Experi- 
ment 2" for this setup is approximately 59 without the 
foreground issues of an autocorrelation. However, this 
would be a much more expensive instrument. 

Instead of using BOSS, we could al so cross-correlat e 
Experiments 1 and 2 with HETDEX (jHill et al.ll200l . 
which will observe 1 million Lya emitters over 200 square 
degrees in the redshift range 1.8 < z < 3.8. In a Az = 
0.007 redshift bin, HETDEX will have an areal density of 
about 57,500 str -1 . This specification gives us an SNR 
for each bin (over full redshift range 2.93 < z < 3.07) of 
5 (22) for Model A and 13 (58) for Model B, which is 
much higher than we can get with BOSS. 

Another option is to just cross-correlate the 
temperature maps from the Planck satellite 
(|The Planck Collaboration! 120061 ) with the same 
photo-quasar map we used in this analysis. However, 
we find that for Model A none of bands have a SNR 
greater than one. We find that the highest SNR we get 
is 0.73 for the CO(2-l) line for Planck's 143 GHz band. 
We assume subtracting the 100 GHz band from all the 
other bands and subtracting the 70 GHz band from the 
100 GHz band. All the other SNRs for the other Planck 
bands are much less for the both lines. Model B is more 
constrainable with SNR=2.0 in the 143 GHz band for 
the CO(2-l) line, but foregrounds will probably degrade 
this signal. Using the BOSS spectro-quasar survey 
would increase / s k y to at most 0.25. Since most of the 
noise errors is due to the temperature maps, Planck x 
BOSS would not do much better than Planck x SDSS 
DR6. Moreover, the use of the CII line instead of CO 
would be more appropriate for Planck high frequencies. 

Current ground-based, high-angular-resolution CMB 
polarimeters offer another promising avenue towards 
measuring this contribution. Consider an SPTPol-like 
survej0 with 6.5 and 4.5 //K-arcmin sensitivity at 90 and 
150 GHz with 2.0 and 1.2 arcmin FWHM and cover- 
ing 500 sq. degrees of a BOSS-likc quasar survey. The 
real SPTPol and BOSS surveys cover opposite hemi- 
spheres, but we perform this exercise as an illustration. 
We assume both bands have a 30% bandwidth. For 
z r-j 3 quasars, the 90 and 150 GHz bands could con- 
strain (Too) for J =3 and 5, respectively. The BOSS 
full survey will have quasar areal densities within the 
relevant redshift bins of about 36,000 and 22,000 str -1 , 
respectively. With these parameters, we forecast cross- 
correlation SNRs for Model A (Model B) of 5 (12) and 

6 gp-ji Team, private communication 
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8 (13). For the extended SPT-3G, the 220 GHz band is 
added, which can constrain the J = 8 line using quasars 
at z ~ 3. The sensitivities for the 90, 150 and 220 GHz 
bands change to 4.2, 2.5, and 4.0 ^K-arcmin and the 
FWHMs change to 1.7, 1.2, and 1.0 arcmin. In the 220 
GHz band, the areal density of quasars in the relevant 
redshift bin is about 15,000 str" 1 . For SPT-3G, we fore- 
cast cross-correlation SNRs for Model A (Model B) of 
15 (31), 24 (31) and 54 (59). HETDEX does not do as 
well as BOSS in this case because it covers such a small 
area. Of course, SPTPol cannot constrain the same lines 
as the 28.8 GHz experiment mentioned earlier. Also, the 
large bands and limiting frequency coverage will make 
the foreground removal difficult so that these numbers 
are optimistic. Also, the high optical depth limit is less 
certain for higher J lines, making the signal more uncer- 
tain for this experiment. Finally, the frequency bands for 
ACPol/SPTPol are wide enough that a full 3D analysis 
is not feasible. 

Other experiments that could also help with this kind 
of search are the P rimordial Inflation Explorer (PIXIE) 
(jKogut et al.l I2011D. and the Mu rchison Widefield Ar- 
ray (MWA) (jLonsdale et al.ll2009h . Although PIXIE is 
a polarization experiment to detect inflationary gravita- 
tional waves, its high frequency resolution over a wide 
frequency range can constrain CII and CO, particularly 
for higher J lines, over large redshift ranges, including 
high redshifts. MWA is currently searching for the 21 
cm HI line from the dark ages and rcionization. As men- 
tioned in Lll, CO x 21 cm can be a powerful probe of 
the high redshift universe. For example, MWA x SPTPol 
could constrain CO J=5 or 7 lines and star formation at 
z 6-7. 



7. CONCLUSIONS 

We have predicted an angular cross-power spectrum 
between CO line emission and quasars and LRGs based 
on ACDM cosmology and the Lll model. We proposed 
searching for the quasar/LRGs cross-correlation to char- 
acterize CO emission in high-redshift galaxies. We have 
also attempted to detect the cross-correlation in WMAP 
and SDSS photo-quasars and LRGs up to z ~ 6. A 
signal was not detectable, mainly due to the large sta- 
tistical errors in the WMAP maps. We were able to set 
upper limits to the brightness temperature of the CO(l- 

0) and CO(2-l) lines, which rule out models much greater 
than our Model B. We also explored the CO(l-0)xCO(2- 

1) cross correlation, another signature of CO emission. 
Although current probes appear to be unable to detect 
CO emission, the potential for future experiments looks 
considerably greater. Current or soon-to-happen ground 
based, high-angular-resolution CMB experiments over- 
lapping with BOSS offer a chance to detect higher J 
lines. In our forecasts for an optimistic model for a fu- 
ture spectrograph to detect CO (1-0) line emission, we 
found a SNR of 58 for a CO(l-0)x(HETDEX Lya emit- 
ter) analysis at z ~ 3 and a SNR of 59 for a more expen- 
sive CO(l-0)xCO(2-l) analysis. Although these numbers 
will likely be decreased due to foreground subtraction, 
this result is still very promising. A future detection of 
CO brightness temperature perturbations will allow us 
to model the CO emission-line galaxies at high redshifts, 
possibly even out to redshifts in the rcionization epoch. 
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